ON TURAN TYPE INEQUALITIES FOR MODIFIED BESSEL FUNCTIONS 

ARPAD BARICZ AND SAMINATHAN PONNUSAMY 



Abstract. In this note our aim is to point out that certain inequalities for modified Bessel func- 
tions of the first and second kind, deduced recently by Laforgia and Natalini, are in fact equiva- 
lent to the corresponding Turan type inequalities for these functions. Moreover, we present some 
new Turan type inequalities for the aforementioned functions and we show that their product is 
decreasing as a function of the order, which has application in the study of stability of radially 
symmetric solutions in a generalized FitzHugh-Nagumo equation in two spatial dimensions. At 
the end of this note a conjecture is posed, which may be of interest for further research. 

1. Some inequalities for modified Bessel functions 

Let us denote with I v and K v the modified Bessel functions of the first and second kind, re- 
spectively. For definitions, recurrence formulas and other properties of modified Bessel functions 
of the first and second kind we refer to the classical book of G.N. Watson [Waj . 

In 2007, motivated by a problem which arises in biophysics, Penfold et al. [PVG] proved 
that the product of the modified Bessel functions of the first and second kind, i.e. u h-> P u {u) = 
I u (u)K u (u), is strictly decreasing on (0,oo) for all v > 0. It is worth mentioning that this result 
for v = n > 1, a positive integer, was verified in 1950 by Phillips and Malin [PMJ. In order 
to shorten the proof due to Penfold et al. [PVG], recently the first author [Ba4] pointed out 
that the Turan type inequalities for modified Bessel functions of the first and second kinds are 
in fact equivalent to some known inequalities for the logarithmic derivatives of the functions in 
the question. For reader's convenience we recall here the historical facts for these Turan type 
inequalities (see [Ba4l lBa7] for more details). More precisely, in view of the recurrence relations 

(1) !„_!(«) = {u/u)I v {u) + I' u (u) 
and 

I v +i(u) = l'„(u) - (v/u)I„(u), 

the Turan type inequality 

(2) (ti)J„+i («) - [I u (u)\ 2 < 
is equivalent to 



2 



(3) ul' v {u)/I v {u) < Vu 2 + v 

where v > — 1 and u > 0. To the best of authors knowledge the Turan type inequality ([2]) for 
v > — 1 was proved first in 1951 by Thiruvenkatachar and Nanjundiah |TN| and later in 1974 
by Amos [Ami p. 243] for v > 0. In 1991 Joshi and Bissu |JB| proved also @ for v > 0, while 
in 1994 Lorch |Loj proved that ([2]) in fact holds for all v > —1/2 and u > 0. Recently, the first 
author [Ba7| reconsidered the proof of Joshi and Bissu [JBJ and pointed out that d2j and ([3|) 
hold true for all v > —1 and the constant zero on the right-hand side of (UJ is the best possible. 
It is worth to mention that in fact the function v i— > I v (u) is log-concave on (— l,oo) for each 
fixed u > 0, as it was pointed out by Baricz [Ba6j . Finally, we note that the inequality ([3|) for 
v > was first proved by Gronwall |Gron| in 1932, motivated by a problem in wave mechanics. 
This inequality in 1950 appeared also in Phillips and Malin's paper |PMj for v = n > 1, a 
positive integer. 
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Similarly, by using the recurrence relations 

(4) K v -i{u) = -( V /u)K v {u) - K' v (u) 
and 

(5) K u+l (u) = -K' u (u) + (u/u)K p (u), 

it is easy to prove (see |Ba41 lBa7] ) that the Turan type inequality 

(6) K v ^{u)K v+l {u) - [K v (u)\ 2 > 
is equivalent to 



(7) uK' v {u)/K u {u) < -^Ju 2 + v 2 , 

where v G R and u > 0. The Turan type inequality © was proved in 1978 by Ismail and 
Muldoon [IMj and recently by the author (see [Ba6, IBa7| ) by using different methods. The 
constant zero on the right-hand side of ([6]) is the best possible. Note that for v > 1/2 the 
Turan type inequality ([6|) appears also on Laforgia and Natalini's paper [LN1]. It is worth 
also to mention here that Ismail and Muldoon [IM] , by using the Nicholson formula concerning 
the product of two modified Bessel functions of different order, proved that [IM] the function 
v i — V K v+a (u) / K v (u) is increasing on R for each fixed u > and a > 0. As Muldoon [Mu] pointed 
out, this implies that v i— > K u (u) is log-convex on R for each fixed u > 0. Recently, by using the 
Holder- Rogers inequality, the first author [Ba6] pointed out that the function v i— > K v (u) is in 
fact strictly log-convex on R for each fixed u > 0. Finally, we note that the inequality ([7]) was 
proved for v = n > 1, a positive integer, by Phillips and Malin |PM] in 1950. 

Recently, motivated by some applications in finite elasticity, Laforgia and Natalini [LN2J 
proved the following inequalities 



and 



[ ' K v ^(u) < u 

where u > and v > in the first inequality, and u > and v E R in the second inequality. 
However, inequality ([8|) is not new, as far as we know it was proved first by Amos [Am] in 
1974. It is important to note here that Laforgia and Natalini in order to deduce ([8]) and ([9]) 
used the Turan type inequalities ([2]) and (|6|). Moreover, inequality ([8]) improves a recently 
deduced inequality by the first author [Ba5j, which is useful in the study of the generalized 
Marcum Q— function applied frequently in radar signal processing. Sec [Ba5, LN2] for more 
details and also [Ball IBa2l lBa3l IB"a6l IBaTl iJBl iLNll lLol ITN] for more details on Turan 
type inequalities. 

Here our aim is to show that the inequalities ([8]) and ([9]) are in fact equivalent to the corre- 
sponding Turan type inequalities ([2]) and ([6]). To see this we just need to rewrite the inequality 
([3]), by using ([1]), in the following form 



I„-i(u) v + Vu 2 + u 2 
I u (u) u 

which is actually equivalent to (jSJ). Thus, the inequalities (J2J), ([3]) and ([5]) are equivalent and 
hold true for all v > — 1 and u > 0. Similarly, in view of ((H), the inequality ([7D is equivalent to 



K„-i(u) > -v + Vu 2 + v 2 
K v (u) u 

which is in fact equivalent to ([9]). Thus, the inequalities ©, ([7|) and ([9]) are equivalent. 
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Moreover, by using the corresponding counterparts of the Turan type inequalities ([2]) and © 
(or inequalities ([3]) and 0) we can obtain in a similar way the corresponding counterparts of 
inequalities ([8]) and Q. More precisely, consider the Turan type inequalities 

(10) ll{u) - I v ^{u)I u+1 {u) < Il(u)l{y + 1) 
and 

(11) K 2 u (u) - K v ^{u)K v+1 (u) > Kl{u)/(l - u). 

Inequality ()10p . which holds for all u > — 1 and u > 0, was proved using different arguments in 
1951 by Thiruvenkatachar and Nanjundiah |TN| for v > —1, in 1991 by Joshi and Bissu [JBJ 
for v > 0, and recently by the first author [Ba3, Ba7j for v > —1. Note that the constant 
l/(v + 1) on the right-hand side of (|1Q|) is the best possible, as it was pointed out recently in 
|Ba7] . The Turan type inequality (jlip . which holds for all v > 1 and u > 0, was proved recently 
in [Ba7j. We note that the constant 1/(1 — v) on the right-hand side of (jlip is the best possible. 
Now, consider the corresponding counterparts of ([3]) and ([7]), namely 



(12) Ul' v (u)/I v {u) > y/u 2 v/(v+l)+V 2 

and 



(13) uK' v {u)/K u {u) > -yjv?u/( v -l) + u 2 . 

The inequality (fT2"j) is equivalent to (fTUj) and holds true for all v > — 1 and u > 0, as it was pointed 
out in |Ba4j . This inequality appears also on Phillips and Malin's paper [PMJ for v = n > 1, a 
positive integer. The inequality (|13|) is valid for all v > 1 and u > 0. For u = n > 1, a positive 
integer, it was proved by Phillips and Malin [PMJ, and for v > 1 real by the first author [Ba7| . 
Note that (|13|) is equivalent to (jlip . which can be verified by using the corresponding recurrence 
relations for the modified Bessel function of the second kind, mentioned above. 

Now, our aim is to show that by using the same argument as in the proof of (JSj) and ([9]) we 
can obtain easily the counterparts of ([8]) and Q. Namely, by using (fTJ we conclude that (fT0|) 
or (fl2|) is equivalent to 



(14) W :< ~ V+ ^ + V \ 

Iu-l{u) —U 

which holds for all u > and all v > 0. Similarly, in view Q we conclude that (jllj) or (|13|) is 
equivalent to 



^-lM — u 

which holds for all v > 1 and all u > 0. We note that inequalities (|14p and (|15p in the above 
forms seem to be new. 

Remark 1. We would like to take the opportunity to point out some minor errors, which we 
found in the papers (Ba4l IBa7| . In [Ba4] the first author claimed that (|13p is equivalent to the 
Turan type inequality [Ba4, Eq. 3.3] 

(16) (y - l)K v ^{u)K v+1 {u) - (2u - 1) [K p (u)} 2 < 



and conjectured that (|16p holds true for all real v > 0. Unfortunately, this claim is not true. 
We note that UfT3\i is actually equivalent to Turan type inequality ([TT]). This can be verified by 
using the recurrence relations (HJ) and ([5]) and by using the fact that K v is decreasing on (0, oo) 
for all v £ M. All the same (|16p is valid for all v E [0, 1] and u > 0, it follows from ([6]), as it 
was pointed out in [Ba7]. Moreover, in [Ba7j the first displayed inequality after Theorem 3.1 
should be 



(17) K 2 u {u) - K u ^(u)K u+1 (u) < -^—Kl(u) 
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just like before Theorem 3.1 in |Ba7j . This is an equivalent form of (|16p and in view of the fact 
that the constant zero is the best possible in ([6]) we conclude that for v > 1 the inequality (fT7|) 
or (|16p does not hold. Moreover, the claim in [Ba7] that (|lip improves (|16p . when v > 1, is not 
true. In view of (|17p . there is no connection between (jlip and (|16p . 

We note that using the recurrence relation K u (u) = K_ u (u), the inequality (fl~6j) can be 
rewritten as 

(18) {y + l)K v _x{u)K v+l (u) - (2u + 1) [K v {u)} 2 > 0. 

In [Ba4] we conjectured that (|18p holds true for all < and it > 0. It is important to point out 
here that in view of the above discussion on (|16p the inequality (|18p holds true for v £ [—1,0], 
but does not holds for i/ < — 1. 

Finally, we note that in the proof of Theorem 3.1 in [Ba7] the explicit formulas for the 
function <f> v : (0, oo) — > R, defined by 

K v -x(u)K v+ x(u) 



4>u{u) = 1 

and its derivative should be written as 

4>v{u) - 



1 


' 4 


u 


.TT 2 J 


4 


POO 




Jo 



uj(t)dt _4_ [°° u(t 2 - u 2 )~f(t)dt 

1 o" 



and consequently 



u 2 +t 2 vr 2 7 (u 2 +i 2 ) 2 
2t 2 7 (t)dt 
(u 2 +t 2 ) 2 ' 

32 f 00 ut 2 j(t)dt 



vr 2 7 (u 2 +t 2 ) 3 ' 

Fortunately, this mistake does not affect the proof of Theorem 3.1 in |Ba7j . the proof of Turan 
type inequalities ([B]) and (|11|) is correct with the above small modifications. 

2. Turan type inequalities for modified Bessel functions 

In this section we are going to present some new Turan type inequalities for modified Bessel 
functions of the first and second kind. As we can see below these inequalities are consequences of 
some more general results on these functions. Before we state and prove the first main result of 
this paper let us recall some basics. By definition, a function / is said to be completely monotonic 
(cm.) on an interval A if / has derivatives of all orders on A which alternate successively in 
sign, that is, 

(-i)"/ (n) («)>o 

for all u E A and for all n > 0. It is known that cm. functions play an eminent role in areas like 
probability theory, numerical analysis, physics, and theory of special functions. An interesting 
exposition of the main results on cm. functions is given in Widder's book [WiJ (see also [MSJ). 
The following results complement the discussion given Section 1. 

Theorem 1. Let I u and K u be the modified Bessel functions of the first and second kind. Then 
the following assertions are true: 

(1) the function v 1— > I \rj}{u) is log-convex on (0, 00) for all u > fixed; 

(2) the function v 1— > K ^{u) is log-concave on (0, 00) for all u > fixed; 

(3) the function v \-± K ^(u)/K ^ + \{u) is decreasing on (0, 00) for all u > fixed; 

(4) the function v 1— > I ^{u)K ^p(u) is log-convex on (0, 00) for all u > fixed; 

(5) the function v 1— >• I ^(u)/K ^{u) is log-convex on (0, 00) for all u > fixed; 

(6) the function v h-> I u (u)/ K u (u) is log-concave on (— 1, 00) for all u > fixed. 
In particular, for all v > and u > 0, the following Turan type inequalities hold 
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I* (u)/K 2 («) < [l^/i^)] [/^faV^^M 



and 



I 2 v+l (u)/Kl +1 { U ) < [I u (u)/K v (u)) [I u+2 {u)/K u+2 (u)] 



Proof. For the proof of (1) - (4) we use the following results for modified Bessel functions: 

(a) the function v 2 i— > I u (u) is cm. on [0, oo) for each u > fixed; 

(b) the function v 2 i— > 1/K u (u) is cm. on [0,oo) for each u > fixed; 

(c) the function v \— > K ^ +a (u) / K ^ +a+n {u) is cm. on [0, oo) for all u > 0, a > and 
n€ {1,2,3,...} fixed; 

(d) the function v 2 \— > I u {u)K u (v) is cm. on (0,oo) for each v > u > fixed. 

Parts (a) and (d) were proved by Hartman and Watson [HW], part (b) is due to Hartman 
[Hal| and part (c) was obtained by Ismail [Isl]. Now, recall that if the function / is cm. and 
the function g is nonnegative with a cm. derivative, then the composite function / o g is also 
cm. (see [Bo] IWi] ). Thus, since v \-± \y/v\' is cm. on (0, oo), from parts (a) - (d) we obtain 
the following results: 

(i) the function v i->- 1 ^{u) is cm. on (0, oo) for each u > fixed; 

(ii) the function v \— > 1/K ^(u) is cm. on (0, oo) for each u > fixed; 

(iii) the function v \— > K ^{u)/K ^ + i(u) is cm. on (0, oo) for all u > fixed; 

(iv) the function v i— )■ I ^(u)K ^{u) is cm. on (0, oo) for each u > fixed. 

On the other hand, each nonnegative cm. function is log-convex (see | Wij ) . and thus parts (i), 
(ii) and (iv) imply parts (1), (2) and (4) of this theorem, while part (3) follows from (iii). 
We note that actually part (3) can be proved directly by using the formula (see [Gros, Is2]) 

Ku-i(V^) _ 4 m_ * _1 



L "^T^' where 



^UK„(y/u) IT 2 J U + t 2 ' '" W J2(t)+y2(t) 

and u, v > 0. Here J v and Y v stand for the Bessel function of the first and second kind, respec- 
tively. More precisely, in view of the Nicholson formula |Wa| 



8 f°° 

J 2 {u) + Y 2 (u) = -o K (2usmht) cosh(2i/£ 
n Jo 



)dt, 

the function v \— > y v (t) is decreasing on (0, oo) and so is the function v i— > 7 v /^ + i(i) for all t > 
fixed. This in turn implies that the function 

K^(u) 4n [°° ly/t+xWdt 



K^ +1 {U) TT 2 Jo U 2 + t 2 

is decreasing on (0, oo) for all u > fixed. 

Finally, we note that part (3) can be obtained also directly by using the fact that the function 
v i y K u+a (u) / K u (u) is increasing on R for each fixed u > and a > (see [IMJ). 

(5) This follows easily from parts (1) and (2) of this theorem. 

(6) Similarly, this follows from the facts (see [Ba6]) that v i— > I u {u) is log-concave on (— l,oo) 
and v i — V K v (u) is log-convex on M for all u > fixed. □ 
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3. MONOTONICITY PROPERTIES OF THE PRODUCT OF MODIFIED BESSEL FUNCTIONS 

Our second main result reads as follows. 

Theorem 2. The function v i-> I v (u)K v {u) is decreasing and (0, oo) for all u > fixed. More- 
over, for all v > 1/2 and u > we have 

(19) 2I u {u)K v {u) < I„_i («)#„_!(«) + («)#„+!(«). 

In particular, the sequence {I n (u)K n (u)} n >i is decreasing and convex and for all u > we have 

(20) Io(u)K (u) > h{u)K x {u) > I 2 (u)K 2 (u) >■■■> I n {u)K n {u) > ■■■ . 

Proof. Recall that the function v 2 i— > I v (u)K v (v) is cm. on (0, oo) for each i> > -u > fixed (see 
[HW]), and consequently the function v i->- J ^{u)K ^{u) is also cm. on (0, oo) for each u > 
fixed (see part (iv) in the proof of Theorem [1]). In particular, the function v i— >■ I /^(u)K ^(u) 
is decreasing on (0, oo) for each u > fixed, and hence the function v \-± I u (u)K v (u) is also 
decreasing and (0, oo) for all u > fixed. Thus, for all u > we have 

h{u)Ki{u) > I 2 {u)K 2 {u) > ■ ■ ■ > I n (u)K n (u) > . 

Now, we are going to prove ([19]) . By using the notation P u (u) = I u (u)K u (u), the inequality 
(|19p can be rewritten as 

(21) 2P u (u) < P v -i(u) + P u+1 (u), 

We note that for the function P v the following recurrence formulas (sec [PVGJ) are valid 

(22) 2uPl(u) = u(P u+1 (u) - P v -x(u)) 
and 

(23) 2vP'l(u) = AuP v {u) - {2u - l)P„_i(u) - (2v + l)P v+l {u). 

On the other hand, owing to Hartman [ Ha2j . we know that the function u \— > uP u (u) is concave 
on (0, oo) for all v > 1/2. Since u \-¥ 2uIi/2{u)Ki/ 2 {u) = 1 — e~ 2 " is concave on (0, oo), we 
conclude that in fact the function u \-t uP u (u) is concave on (0, oo) for all v > 1/2. Consequently, 
for all u > and v > 1/2 we have 

uP'M < -2P' v (u). 
Now, combining this with (|22|) and (|23|) we obtain for all u > and u > 1/2 

2v \2P v (u) - {P u -i{u) + Ph-i(«))] < 0, 
which is equivalent to (|2ip . Finally, by using (|19j) for y = lwe obtain 

Pi(u)-P («) <P 2 («)-Pi(«) 

and hence by using the fact that the sequence {I n (u)K n (u)} n >i is decreasing for each u > 0, it 
follows that the sequence {/ n (u)i ; C ri (?/)} n >o is also decreasing for each u > 0, i.e. (f20j) is valid. 
With this the proof is complete. □ 

Comments and concluding remarks 

1. In j DHK] IHDK} HDKP (see also [He]), the authors study the existence, stability and 
interaction of localized structures in a one-dimensional generalized FitzHugh-Nagumo 
type model. Recently, van Heijster and Sandstede [HS] started to analyze the existence 
and stability of radially symmetric solutions in the planar variant of this model. The 
product of modified Bessel functions I u {u)K u {u) discussed in this note arises naturally 
in their stability analysis, and the monotonicity condition (|20p is important to conclude 
(in)stability of these radially symmetric solutions. 
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After we have finished the first draft of this paper, van Heijster informed us that the 
monotonicity of v i— )■ I u (u)K u (u) is in fact an immediate consequence of the integral 
formula [MPS! p. 98] 

/•oo 

I u (u)K u (u)= / I 2v (2usmht)e- 2ucosht dt 
J o 

and the monotonicity of v i— )■ I v {u). More precisely, due to Cochran [Co] we know that 
v I—)- I„(ii) is strictly decreasing on [0, oo) for all u > fixed. By using this result and the 
above integral formula we conclude that in fact v i— > I u (u)K u (u) is strictly decreasing on 
[0, oo) for all u > fixed. Thanks are due to van Heijster for the above information. 

2. Let us recall that the function u h- > P u (u) = I u (u)K u (u) is strictly decreasing on (0,oo) 
for all v > —1. For, v = n > 1, a positive integer, this was proved in 1950 by Phillips 
and Malin [PM]. for v > real in 2007 by Penfold et al. [PVGj . and for v > -1/2 and 
i> > — 1 recently by the first author |Ba4[ IBa7j . We note that by using ([22]) and the 
monotonicity of z/ i— > P u (u), proved in Theorem[2J we obtain immediately that u \— > P u {u) 
is strictly decreasing on (0, oo) for all u > 1. 

Now, let us consider the function u i— > 2uP u (u). Hartman and Watson [HW] proved 
that this function is a continuous cumulative distribution function on (0, oo) for all 
v > 1/2. This function actually maps (0,oo) into (0, 1). Moreover, later Hartman [Ha2j 
proved that u h-> 2uP u (u) is concave on (0, oo) for all v > 1/2, and as we pointed out in 
the proof of Theorem [21 this result is valid also for v = 1/2. We would like to point out 
here that actually the monotonicity of u i— > P u (u) for v > 1/2 is almost trivial by using 
the aforementioned concavity result of Hartman. To see this let us recall a particular 
form of the Pinelis version of the monotone form of l'Hospital's rule (see [Pi]): 

Let f,g : (a,b) C R — > R be differentiable functions on (a,b) with g'(u) ^ for 
all u G (a, b). Furthermore, suppose thai lim n _> a f{u) — lim u _^ a g(u) = and f'/g' is 
decreasing on (a, b). Then the ratio f/g is decreasing too on (a, 6). 

Appealing to this result, since lim u _>o 2uP u (u) = lim u ^o u/v = 0, for v > 1/2, to prove 
the monotonicity of u t-t P u (u) = uP v {u)/u it is enough to show that u i— > [uP v (u)]' /u' 
is decreasing, which is clearly true because u \— > uP u (u) is concave. 

3. Finally, we note that the product P v is useful also in other problems of applied mathe- 
matics. For example, in 1986 Cantrell [Ca] derived tight upper bounds for the function 
u i y uI n+ i(u)K m+ i(u) , in order to obtain suitable truncation and transient errors in the 
computation of the generalized Marcum Q— function. 

Motivated by Theorems [T] and [2] we conjecture the following. 

Conjecture. The function v i— > I u (u)K v (u) is log-convex on (— l,oo) for all u > fixed. 
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